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Abstract 



There is a natural way to study the long distance interactions of gauge 
theories in the electric (momentum) representation. Here, the main ideas are 
presented for the Abelian and Yang-Mills gauge theories emphasizing on the 
structure and the advantages of this approach. 
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I. INTRODUCTION 



For simple 1 + 1 dimensional theories the information provided from the existing symme- 
tries is enough to determine the dynamics of the theory (wave functional, correlators etc.). 
Such examples are the conformal theories |l[ as well as the pure gravitation []2],|3| and Yang- 
Mills [|]|| , which in 1 + 1 dimensions are easily solved just by considering the constraints 
the (conformal, general coordinate or gauge) symmetries impose on the wave functional. In 
these low dimensional field theories, except in the case we incorporate additional fields, the 
Hamiltonian is trivial (e.g. has only kinetic term) and it is not difficult to find its eigenvalues. 
Hence, the information of the constraints applied on the Hamiltonian through the properties 
they impose on the wave functional is enough to determine completely the dynamics. 

In higher dimensions the Hamiltonian is more complicated PHJlH. In order to solve the 
theory apart from solving the constraints corresponding to the symmetries of the theory 
you need to perform non-trivial calculations for determining the energy eigenvalues and the 
corresponding wave functionals. This does not diminish the contribution of the information 
provided from the constraints. In fact we shall see that in the electric representation the long 
range interactions are determined mainly by them through the implementation of particular 
constraints they impose on the field configurations. These, by their turn, incorporated in 
the Hamiltonian eigenvalue equation will give the particular interactions dictated from the 
theory. 

The advantage of working in the electric (E) representation |5|,p^-[l9[| (rather than the 
vector potential (A) representation [|,^-|23|) is that by separating the gauge degrees of 
freedom as the ones which rotate the E vector the constraint becomes a multiplicative (non- 
differential) equation. Fixing the gauge within the above procedure by fixing the direction of 
E, a part of the gauge degrees of freedom are removed. The remaining ones, as we shall see 
for the case of the Yang-Mills and Abelian theories, are the rotations around the direction of 
the electric vector, E. Integrating these angles out constrains the field configurations via the 
generation of delta functionals. As a result the quantum constraints are reduced to classical 
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ones, which are imposed as the arguments of these delta functionals. Still the procedure here 
is in some aspects simpler than the classical case: it is not necessary to solve the equations 
constraining the field configurations, but rather they can be used as definition equations for 
the divergence of the field, diE 1 , and can be substituted in the Hamiltonian when acting 
on the wave functional. In this way the terms in the Hamiltonian eigenvalue equation are 
reduced significantly 

To understand the above procedure, first we shall rigorously study the Abelian theory 
in the E representation of the field within the Schrodinger representation including static 
sources. There, the usual results of the vacuum energy and the Coulomb potential between 
the sources are derived [ 24| . In addition the Fourier transformation of the E representation 



vacuum wave functional is performed, and the corresponding one in the A representation is 
derived where the sources appear naturally to be "dressed" p5| , p6| . Then we turn to the 



Yang-Mills theory with gauge group SU(N), where by proceeding in the same way as in the 
Abelian theory the Gauss law constraint is solved, which imposes Abelian like constraints 
to the N — 1 Cartan components of the electric field. The implications for the generation 
of the Coulomb and confining interactions are presented without solving the Hamiltonian 
eigenvalue problem JHj]. 

For an introductory study of Schrodinger representation see []27H3lH . 



II. ABELIAN THEORY 



As a first step in order to understand the E representation we should solve the Abelian 



theory, where we incorporate only static sources rather than fermionic fields ||26|| . In the 
E representation the field variables are E l (x), while A l (x) = i5/5E l (x). For the wave 
functional of the theory, ^[E], the Gauss law with two sources of opposite unite charges is 

diE\x)^[E] = (5(x - x ) - 5(x - . (1) 



2 



As the operator diE l (x) acts multiplicatively on the wave functional, equation (|]) can only 
be solved by imposing on ty[E] to vanish everywhere else accept when the field configuration, 
E, satisfies 

diE l (x) = 5(x — x ) — 5(x — xi) , (2) 

which is the classical Gauss law of the theory. Hence the wave functional can be taken to 
have the form ty[E] = Y[ x 8\diE i (x) — 5(x — xo) + S(x — xi)^J§[E}, where $[J5] is a gauge 
invariant functional. This is the proper way to impose such a constraint in quantum field 
theory, where E should take all possible configurations. 

There is another way to approach the Gauss law. Equation ([]]) determines that under a 
gauge transformation ty[E] should transform in the following way 

V[E] -» ^[E u ] = e^ ixE ^ u ' lu U{x l )m[E}U-\x ) = 

e -i J d 3 xE i (x)a i e(x)-ie{xo)+ie(x 1 )^j^ ^ 

where U{6) = expi9 is a U{\) gauge group element. Note that E u = UEU~ l = E, 
i.e. in the Abelian theory the electric vector, E, is gauge invariant as it does not change 
under a 9 rotation. So the wave functional $?[E] should not get transformed under U. The 
phase factor in the transformation is produced from the phase factor exp — i J E l A l when 
the Fourier transformation of the gauge invariant functional *&[A] is taken (for more details 
see H). This phase factor should be combined with the group elements U(x\) and U~ 1 (xq) 
and cancel out, which happens when the condition (|2]) is imposed on the field configurations. 
Hence, the symmetry of the theory straightforwardly imposes the above constraint. Also, 
from (0) we see that the solution of (|1|) can be expressed as the group integration over [/(!):[] 



x This integration over the gauge degrees of freedom is similar to the procedure followed in [25|. 
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V[E] = J Vue-I d3xEi9iUU 'uixJu-^Xo^iE] = 
J V9e~ i f d3xEi ( x ) d i 6 ( x )- i6 ( x o)+ iS { x i)^jr;^ = 

Y[S(diE l {x) - 5 {3 \x - x ) + S^ix - . (4) 

X 

That is the integration over the gauge degree of freedom 9 gives the delta functional con- 
dition. The functional integration T>9 is along the range [0, 27r], but it can be extended 
to [— oo,+oo] by assuming that the charges of the theory are quantized, which makes the 
exponential periodic in 2ir, as it should be 0. Though, there is an overall factor inserted 
from this extension, which makes the delta functional approach one rather than infinity 
when its argument goes to zero. For this reason we shall take the action of the functional 
differentiations on this delta functional to vanish. 

§[E] will be determined from the following Hamiltonian eigenvalue equation 



HV[E] = - f d 3 x[g 2 E i {x)E i {x) + \B i (x)B i (xj)^[E] = £V[E] (5) 
2 J v Q 



where 



B\x) = &%A\x) = ^%J^ ■ (6) 



The Hamiltonian in the momentum representation of the spatial coordinates takes the form 

H — — I -p) - -^e ijk e ffk [ d 3 p(27t) 3 p lPl -4rT — S — - . (7) 

2 J (2vr) 3 KF > v y> 2g 2 J yK ' FjFj 5E k (p) 5E k {-p) KJ 

Here the Fourier transformations are defined as 

E'(*) = J^'~E'(P) -d = / d 3 p(2x)V ^_J_ . (8) 

2 A requirement, which would be enough to impose the constraints by itself. 
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We can proceed in analogy with the A representation where the vacuum wave functional is 
exponentially quadratic in the magnetic field by defining a new "magnetic" field composed 
from the electric (multiplicative now) variable, E, as 

B\x) = ie ijk d j E k (x) . (9) 

Hence an ansatz for the unknown functional (in the momentum space) is 

m =exp / (^(^(p^HO = 

exp - J TjSgM*) (k 2 E\k)E\-k) - k l E l (k)k j E j (-k)) (10) 



where we have used the relation e^ k e^ k = 5^5 kk — 5^ k 5 k K The function h(p) is unknown. 
The full wave functional, ^[-E], includes the delta functional mentioned before, while its 
condition can be now applied on the exponential to simplify it by substituting the second 
term in the exponent with the delta functions denoting the positions of the sources. But 
on this exponential will act functional differentiations so they will give zeros when they see 
these delta functions; so we can neglect them without any loss of information. As a check 
we could keep in the exponential all the terms as in (|1(]) deriving at the end exactly the 
same result obtained from the simplified functional Q. From now on we shall disregard these 
terms as we do not want the delta functionals to act on $>[E], but rather on the eigenvalues 
of the operators we are interested into. 

Applying the Hamiltonian operator on ty[E] we obtain 

2^ / (2^ 4p4 h2 (p)(p 2E Kp) E X~P) -P iE W EJ (-P))]nE] (11) 



3 A somewhat expected result as k l E l (k) is longitudinal and the Hamiltonian should be "trans- 
parent" to these terms. 
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Obviously the third term in the RHS can be made to cancel the second by choosing, h(p) = 
g 2 /(2\p\ 3 ). But there is still dependence of the eigenvalue, S, on the field from the fourth 
term! At this point the Gauss law comes to rescue with the delta functional constraints, 
which impose the condition (0). Hence, we derive 

' 2 ' d'p 1 



£ = 5 (3)( ) / d 3 p M + y/ j^^p'E^E^-p) 



<$< 8 >(0) J d 3 p \p\ + jj d 3 xd 3 y J i e -*(3H0 <),!■;■ (,■)<), I-yijj) 



5 (3) (0) / d 3 p \p\ + 



g ! f d 3 xd 3 y — -^—(5 {3) (x-x )-S {3) (x-x 1 ))(S {3) (y-x )-S^(y-x 1 

2 J 47T \X — V \ X ' V 



e = s v i0) J d s pM _f_L_L_ + vii (12) 

where V®^ is the self energy of the sources. This is the well known result. The energy of the 
vacuum as well as the attractive Coulomb potential between the two sources is also easily 
derived in the A representation. 

We see that \fr[E] is exponentially quadratic in the field so it is possible to perform its 
Fourier transformation. For this we rewrite the delta functional as 

HsfdiE^x) -5 {3) {x-x ) + 5 {3) {x-x 1 )) = 

X 

J VOexpi J J^Lfa&ik) - e-**° + e-**)9{-k) . (13) 
Taking the Fourier transformation of ^[-E], we derive 



9[A] = J VEe W {-z J ^E\k)A\-k)}m[E] 
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f VEVOexp [ -^-\ - i£-LE i (k)E i (-k) - E\k){iA\-k) + fc^(-fc)) 

J J (271") 2 rC 



(27T) 



f(e-**' 3b - e~ ik - 3l )9{-k) 



[V6exp — \ f [|jfe|4*(jfe)4*(_jfe) -2i\k\k i A i (k)9(-k) + \k\ 3 9(k)9(-k) 
i2g 2 {e- i% - 3 ° - e~ itSl )9{-k) 

9[A] = exp { - ^5 / ^y ^y - zA(z ) + iA^)} , (14) 

where B l (x) is the magnetic field constructed out of the multiplicative now A l (x) (see |^4}|), 
while 9jA(x) is the longitudinal component of A 1 (re) with 

A(,)= iffytm. (15) 

The magnetic term B l (x)B l (y) can be also written with respect to only transverse compo- 
nents of A. The exponential of A(x) as presented in the wave functional, can be combined 
with the fermionic sources denoted by operators multiplied on *&[A] making them "dressed", 
that is gauge invariant! The dressing of the sources is necessary to give the proper asymp- 
totic behavior to the fermions |[26|| . However, performing the dressing with a Wilson line of 
the full A vector rather than only its longitudinal component, would result to infinite terms 
in the potential as explained in ||33|| . That the delta functional evokes only the longitudinal 
component of A can be seen by Fourier transforming the delta functional Y[ x ^{diE 1 ) alone 

J VEV9exp-i J -^-E^k^A^-k) - ik l 9(-k)) = J V0S(A*(x) - di0(xj) (16) 

which picks up only the longitudinal components of the vector potential for the dressing of 
the sources. 

In addition, the simplification of the functional in the E representation due to the delta 
functionals is not possible in the A representation as there isn't such a condition there 



(the Gauss law is implied via the gauge invariance and A(x) terms for the sources). This 
simplification may be of small importance for the Abelian theory, but in the case of the 
Yang-Mills theory, where the terms needed to consider in the wave functional are many, a 
reduction of this form may be of significance. 



III. YANG-MILLS THEORY 

Here, we shall study the more intriguing problem of Yang-Mills theory in the E repre- 
sentation [|l2] , |l~T| , [3~8j| . The classical Hamiltonian is given by 



r a 2 1 

H = J d 3 x^E ia (x)E ia (x) + —B ia (x)B ia (x)) , (17) 

where B ia (x) = e^ k {djA ka {x) + \j2f abc A^ b {x) A kc (x)), and it is accompanied by the Gauss 
law constraint 

diE ia {x) + if abc A ib {x)E ic {x) = , (18) 

where the index a runs over the N 2 — 1 components of the su(N) algebra. For quantizing 
in the electric Schrodinger representation the commutation relation 

[A ia {x), E jb (y)] = i5 ij 5 ab 5^(x - y) (19) 

is materialized by taking E m diagonalized and A ia a differential operator, A ia (x) = i SE f a ^ ■ 
The wave functional, ^[-E], has to satisfy the Hamiltonian eigenvalue equation H^/[E] = 
£ty[E] and the Gauss law constraint 

(diE^ix) - ^E^jJ^y) ^ = nE]T a 5^(x - x Q ) - n[E}5^(x - x x ) . (20) 

where a source is placed at x and an antisource at X\. Hence, the Gauss law enforces the 
following transformation property on ^[i?]: 

m[E u ] = e^f^'^Uix^iEp-^xo) , (21) 

similar to the Abelian case, but now U G SU(N) is a matrix. The solution of ( PUD can be 
given in terms of the following group integration 
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V[E) = J Vue-^ r I Ei9iUU "uixt^ixo) $[E) (22) 

The undetermined functional, <&[E], is invariant under gauge transformations and will be 
calculated from the Hamiltonian eigenvalue equation. In contrast to the A representation, 
where gauge invariant objects made out of A are non-local, for example the trace of the 
Wilson loop, in the E representation the gauge invariant objects are more elementary due 
to the vector character of the field. They are traces of products of the field E. 

We can separate the gauge degrees of freedom from the dynamical variables ||. As the 
gauge transformations rotate the E vector in the su(N) space we can consider one fixed 
direction, named K l ; then the general vector E % is obtained from it by an SU(N) rotation, 
g, as E l (x) = g(x)K t (x)g~ 1 (x), where g transforms as g(x) — ■> g u (x) = U(x)g(x). Hence, 
the SU (N) gauge symmetry of the electric field is clearly separated. K l can be used as the 
dynamical variable while g can be fixed as its value should not have any physical relevancy. 
K ia has 3(iV 2 — 1) degrees of freedom and g has iV 2 — 1. So the reparameterization of the 
electric field in terms of K and g needs additional constraints in order to become one-to- 
one. To define the appropriate constraints we need particular decomposition properties of 
SU (N) , which have some interest by themselves. 

A general SU(N) element, g, can be taken to satisfy 

J a L g = L ab ^- b g = -T a g , J a R g = R ab J^g = gT a , (23) 

for a special reparameterization, g(x), with respect to the N 2 — 1 angles, Xa, and for some 
invertible matrices L and R ||. From the N 2 — 1 matrix hermitian generators, T a , of the 
SU(N) group, satisfying tiT a T b = c5 ab , it is convenient to let the Cartan elements to be 
the N — 1 first ones, i.e. {T K ) a p = f K (a)S a ^ for k = 1...N — 1. Along with this matrix 
representation of the generators the differential generators or J R can be organized 
similarly by diagonalizing the upper (N — 1) x (N — 1) block of L and R. The first N — 1 
of the jT's can be taken to be the differentiation with respect to a single angle, named 4> K 
for the "left" generators and K for the "right" ones. Then the first iV — 1 of the relations 
(g3|) become 
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^ia = ~i^9 = -T K g , J K R g = -i~g = gT\ «=1..JV-1, (24) 

which can be solved to find the <fi and dependence of g, as 

<?(x) = h{<f>)~g(9)Ht) (25) 

where # are the remaining (N — l) 2 angles of x, and the h elements belong to the Cartan 
subgroup, H, of G = SU(N), while g G G/H. The range of the angles is <fi : [0, 27r], 
4> : [0,2tt] and 9 : [0,tt]. 

With the above diagonalization of the SU(N) group it is possible to perform the group 
integration in relation fl22"|). The functional turns out to have only exponential dependence 
on 0, so the T>(j) integration produces a delta functional. (The angles are similar to the 
9 angle we met in the Abelian case. For the algebraic components of K l which belong in 
the Cartan subalgebra the group elements h(<f>), in relation gK l g~ l , commute with it and 
cancel out! These N — 1 components, K tK , turn out to have Abelian-like behavior). The 
generated delta functionals make the T>9 and T><p integrations easy to perform, resulting 



finally in yjjrg 



q?[E] = e-^I^-'^g^P'g-'ixo) 



II U6(d i K iK (x)-r(p)5W(x-x ) + r(p)6®(x-x 1 )) $[E] , (26) 



N-l 

X 

K=l X 

where P p = diag(0...0, 1, 0...0), with 1 in the p-th place, while its dimensionality depends 
on the representation of SU(N) we have chosen. The first two terms of (|26"D give the 
transformation properties of the wave functional required from the Gauss law. In the case 
where there are no sources, the constraints in the delta functionals are diK lK = for k = 
1...N — 1. It is natural then, to extend the constraints to diK m = for a = 1...N 2 — 1 
so that spatial symmetry is restored and the decomposition E % = gK % g~ l becomes one-to- 
one ||34j| . This symmetry is obviously broken when the sources are present. We shall keep 
the constraints for the rest of the components in mind without inserting them in a delta 
functional form. 
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A similar argument to the Abelian one for the generation of the constraints by conditions 
imposed from the phase factor produced from a gauge transformation cannot be applied here 
as a simple calculation shows that the phases produced from (p6|) cancel identically. This is 
due to the existence of expifl = exp — -tr / E l d{gg~ x in the wave functional, which creates 
after a gauge transformation the proper phase factors for the cancelation. Note that Q is not 
an analytic functional of E, as it depends also on g, and hence you cannot write Q[E}. Still 
the demand of periodicity of the angles in Q generates the above mentioned constraints 
on d l K iK . 

In order to determine <&[E] the Hamiltonian has to be applied on ty[E] and its eigenvalues 
should be sought. For this it will be necessary to fix the symmetry of the rotations with 
respect to g G SU(N), as these variables will cause extra divergences when the functional 
differentiations are taken |L^,[L4]]. These divergences are not relevant for the dynamics of 
the theory because the physical variables have to be independent from the gauge rotations. 
Hence, we take g — 1, and the wave functional becomes 

AT-l 

V[K\=Y,P p II Yl^K^tx) - t(p)5^\x - x ) + r(p)5^\x - Xl )) $[K] . (27) 

P K = l X 

The presence of the sources is still denoted with the delta functionals. The unknown part is 
the gauge invariant functional 3>[-K], which will be determined from the Schrodinger equation. 
Note the similarity of this functional with the one in the Abelian case (|j). With the electric 
field fixed towards the K l direction the Schrodinger equation becomes 

r n 2 1 
H[K\V[K\ = / d 3 x(— K ia (x)K m (x) + —B ia (x)B ia (x))^[K] = S^[K] (28) 



with 



Si & & 

B™(x) = ^ k (dj 6Kka ( x) + 2 fabC 5K^(x) 5K^{x) ' (29) 

Here we shall not attempt to solve this equation. Though, it would be tempting to try 
an ansatz as in the Abelian case to be the exponential of / d 3 xd 3 yh(x,y)B m (x)B ia (y) for 
B being the modified magnetic field constructed out of iK ta (x) in the place of A m (x). 
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Merely, we shall be interested in determining the conditions for deriving the Coulomb and 
the confining potentials with the help of the delta functionals generated from the Gauss law 
or equivalently implied from the gauge symmetry of the theory. 
Let us take the following functional 

*[*!= «P {A„ / ^ m **<*)**<-*) + A. / ^SW^(m"(-k) + ...} - 



= expW eff [K] , (30) 

where W e ff is a proposed effective functional obtained after the renormalization of the wave 
functional is performed (see |7, 32, 35, 36, 17]). There are more terms necessary to fulfill the 
Schrodinger equation but we shall consider only these two as they will produce the potentials 
we are interested into. 

First note that the constant A is dimensionless, so very naturally this term could exist 
in the initial wave functional before renormalization. This term is the same as the Abelian 
one, and its existence in ty[K] can be proven also perturbatively, having in the lowest order 
Ao = —g 2 /2, [0. In the second term, the constant Ai is dimensionful (Ai ~ and this 

should be a result of the regularization of the theory, which generates dynamically mass, 
m, as there is not any dimensionful constant in the initial theory. Hence, this constant, 
Ai, can be written in terms of m. Its calculation would need to consider higher orders in 
perturbation theory. 

For the following we shall also consider only the first term of the potential of the Hamil- 
tonian (p8|) 

^l Wpi sih/iKk=i) (31) 

which as we have seen in the Abelian case gives the Coulomb potential. Due to its structure 
it generates longitudinal components of the electric field, which combined with the delta 
functional conditions produces the interactions. When it is applied on the first term it will 
give in the same way as before the desired Coulomb interaction between the static sources 
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present also in the Yang-Mills theory. In fact, due to the matrix form of ^f[K] we have to 
take the trace of the expectation value of the Hamiltonian. Eigenstates of the Hamiltonian 
would be each term of the sum yj p in (p7|), while the state we are interested into is their 



superposition, which has the desired expectation value (explained in detail in |0|). In 
addition a summation has to be performed over the N + 1 directions the Cartan subalgebra 
can take in su(N), which is a remnant of the group integration (^2|). Hence, the Coulomb 
potential, Vcouh coming from this term is given from 

where C^iV) is the quadratic Casimir invariant of SU(N) and V, ; 1 is the self energy of the 
sources due to the Coulomb interaction. 

The application of ([TI|) on the second term of W e ff will give, up to an additive constant, 
2A ? f j*„ 1 



J d z V - Vl K^[p) V3 K^{- V ) 



g 2 J p 



^$ J d 3 xd*y J ^^e-^^d.K^djK^y) (33) 



? J "J (2vr)3p' 

We can easily calculate that 

d3 P 1 ^iv-(x-v) 1 



sin e cos e n 



where e is going to zero. Hence, the term we are interested in becomes after averaging, 
tracing and multiplying by iV + 1 

Vconf = ^C 2 (N)\ Xl - x \ + V*i , (35) 

where V^ ; 1 is the self energy of the sources due to the confining interaction. It is worth 
noticing that this self energy (which is an infinite term) is produced from the evaluation 
of the integral in ( p4|) rather than the coincidence of the delta function arguments, which 
produces the self energies for the Coulomb interaction. This result reveals the dual character 
of the confining potential ~ \x — y\ with respect to the Coulomb one ~ l/\x — y\. For the 
second the infinities come from the identification of the points x and y while in the first case 
they come from the region around p — > or equivalently x — > oo. 
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IV. DISCUSSION 



We have seen that Yang-Mills theory in the E representation, with the expression of the 
Gauss law with classical-like constraints obtains N — 1 Abelian-like variables: the Cartan 
components of K l . Moreover, the importance of the longitudinal Cartan components of the 
electric field is apparent for the propagation of the interaction between the sources. Their 
divergence, diK lK , is constrained with delta functionals in the same way as in the Abelian 
theory, resulting to their decouple from the Hamiltonian eigenvalue problem. This might 
bring the problem to a simpler form. Moreover, the classical constraints are just used as 
definition equations for d{K lK , so you do not need to solve them and deal with classical 
configurations like monopoles, instantons and so on. 

From this treatment we see that the quantum theory encloses in some sence a classical 
one. This could be connected with the success of the large N study of gauge theories where 
the dominant classical behavior gives the basic characteristics of the theory. Even more there 
have been evidences for some theories that their results do not change much if you take N 
to be quite small, signifying the dominance of the classical character in the full theory |37| . 

We have seen how important is for the calculation of the confining potential to perform 
first the renormalization of the theory, where a dimensionful scale is generated. This will be 
used to change the conditions from the dimensionless factor of the Coulomb term, l/\x — y\, 
to the confining one, \x — y\, which needs to have a dimensionful factor (~ [£] -2 ) in order 



to be in dimensional agreement with l/\x — y\, as seen in (j35|) . 

In the A representation the dressing of the sources in the Yang-Mills theory is not a 



straightforward problem |26| . In the E representation the insertion of the sources is denoted 
simply with their presence in the delta functionals. As we have seen for the Abelian case 
the dressing in the A representation can be obtained with a Fourier transformation of the 
functional ty[E]. But this is hard to perform in the Yang-Mills case where the logarithm of 
$>[K] is of higher order than quadratic in K and also undetermined! So this is an additional 
reason for working in the E representation for the study of fermion sources. 



14 



Of course there are disadvantages in the study of Yang-Mills theory in the E represen- 
tation. In contrast to the kinetic term in the Hamiltonian, which is quadratic in E, the 
potential is up to quartic in A so there are present quartic functional derivatives! This in 
the first sight complicates extremely the problem. Though, an approximate solution of the 
Schrodinger equation focusing on the zero propagating modes is possible with the consid- 
eration only of quadratic terms ||17| . Its presence is enough to guaranty the existence of 
the main features of the theory like mass generation, Coulomb and confining potentials. 
Though, there is need of farther study of the dynamics of the theory and its behavior to 
achieve better quantitative results. 
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